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Abstract. I review the linearized hydrodynamical treatment of a fast parton traversing a perturbative 
quark-gluon plasma. Using numerical solutions for the medium's response to the fast parton, I obtain the 
medium's distribution function which is then used in a Cooper-Frye freeze-out prescription to obtain an 
azimuthal particle spectrum. Two different freeze-out scenarios are considered which yield significantly 
different results. I conclude that any meaningful comparison of azimuthal hadron correlation functions to 
RHIC data requires implementing a realistic freeze-out scenario in an expanding medium. 

PACS. 12.38.Mh Quark-gluon plasma 



1 Introduction 



An exciting result from heavy-ion experiments performed 
at the Relativistic Heavy-Ion Collider (RHIC) is the cre- 
ation of a state of deconfined quarks and gluons known as 

■ the quark-gluon plasma (QGP) (see, for example, [1]). Fast 
] partons (partons which propagate through the medium 
. with a velocity approaching the speed of light) are cre- 
' ated by hard transverse scattering in the early moments 
, of a heavy-ion collision and are a useful probe of QGP 

■ physics. As a fast parton propagates through the medium 
I it interacts with the surrounding plasma, depositing en- 
i ergy and momentum through the process of jet quenching 
' (see, e. g., [2l[3]). A natural question arises: how does the 

energy and momentum deposited by a fast parton affect 
the bulk behavior of an evolving QGP? 

This question has become particularly interesting in 
light of the results of experimental measurements at RHIC 
[HISlinilT] of hadron correlation functions which suggest the 
fast parton may produce a propagating Mach cone in the 
medium. The hadron correlation function is a measure- 
ment of the azimuthal dependence of hadrons associated 
with a high energy trigger hadron. For instance, two par- 
tons may experience hard back to back transverse scat- 
tering near the surface of the medium. One parton will 
rapidly escape the medium and decay into a jet of hadrons, 
one of which is measured and serves as the near side trig- 
ger. The second parton propagates through the medium in 
the direction of the away side, depositing energy and mo- 
mentum. The hadron correlation functions show a double 
peak structure in the away side distribution which is con- 
sistent with the formation of a Mach cone-shaped emission 
pattern. 



Several authors have done a theoretical investigation 
into how the energy and momentum deposited by the fast 
parton affects the bulk behavior of an evolving QGP (see, 
for example, [8l[9l[T0l|TT|[T2l[T3l[TWT5TO ) ■ In light 

of the strong evidence [201121] that the QGP produced at 
RHIC obeys the hydrodynamic assumption of local ther- 
mal equilibrium, the common theoretical approach has 
been to treat the fast parton as a source of energy and 
momentum coupled to the hydrodynamic equations of the 
medium. This type of analysis requires deriving a form 
for the hydrodynamic source term generated by the fast 
parton. In a recent work, Betz et al. j22j made a direct 
comparison of two hydrodynamic source terms, one de- 
rived in the context of perturbative quantum chromo- 
dynamics (pQCD) [18], and the other calculated in the 
strongly coupled Ads/CFT correspondence [15]. In their 
comparison, the authors solved for the hydrodynamic evo- 
lution of an infinitely large, static background medium in 
the presence of the source terms. They then obtained an 
azimuthal away side associated hadron correlation spec- 
trum after performing an isochronous Cooper-Frye (CF) 
freeze-out. The authors found propagating Mach cones in 
the hydrodynamic evolution resulting from both source 
terms; however, in the final hadron correlation spectrum 
the conical Mach structure was washed out by the CF 
freeze-out scenario, in disagreement with the experimen- 
tal data. An exception was found in flow generated by the 
non-equilibrium Neck zone in the Ads/CFT case, which 
allows for a conical structure in the final hadron correla- 
tion spectrum; however, this contribution does not obey 
Mach's law. 

While the work discussed above allows for an interest- 
ing comparison of two different source terms in a particu- 
lar setting, one has to question whether an isochronous CF 
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freeze-out provides an accurate description of the hadroniza- 
tion that occurs in heavy-ion colhsions. An isochronous, 
or fixed time, freeze-out assumes that all of the matter 
hadronizes at the same time. The freeze-out prescription 
requires integrating the distribution function over the freeze- 
out hypersurface, which, for an isochronous freeze-out, is 
just the entire volume of matter. This prescription is par- 
ticularly questionable for the case of a fast parton propa- 
gating through the QGP at RHIC, where it is likely that 
a Mach shock would push matter out of the medium. This 
matter would then hadronize before the rest of the matter 
in the medium. The effect of integrating the distribution 
function over the entire volume of matter at a fixed time, 
is that any cylindrically symmetric, or conical, contribu- 
tions tend to be washed out. In what follows, I will exam- 
ine the azimuthal associated hadron correlation spectrum 
resulting from the hydrodynamic evolution generated in 
the presence of the pQCD source term derived in [l^. It 
will be shown that the final spectrum is sensitive to the 
specific freeze-out assumptions employed. I conclude that 
any meaningful comparison with experimental azimuthal 
hadron correlation functions measured at RHIC requires 
implementing a realistic freeze-out scenario in an expand- 
ing medium. 



2 Hydrodynamic Evolution 

I now consider the hydrodynamical equations for a system 
in the presence of a source term, , which are 



r 



(1) 



where T'*" is the energy-momentum tensor of the system. 
The equations can be linearized provided that the energy 
and momentum density deposited by the fast parton is 
small compared to the equilibrium energy density of the 
medium. In this case one can write T'"' = T^^ + ST'^" , 
where ST^'^ is the perturbation of the energy-momentum 
tensor resulting from the source in an otherwise static 
medium. Eq. ^ now takes the form 



where ^^,T^'' = and (5T^'' is given by [8] 

EE 5e, = g, 



(2) 



(3) 



In ([3]), 5e and g are interpreted as the energy density and 
momentum density, respectively, deposited by the source 
parton, which interacts with the medium via J''. Also, 
Cs denotes the speed of sound, Fg 



At] 



in 



3(eo-l-po) 3sT 

the sound attenuation length, eo and po are the unper- 
turbed energy density and pressure, respectively, and rj is 
the shear viscosity of the medium. 

5e and g must now be solved for in terms of J'^ . It is 
easiest to work in momentum space by taking the Fourier 



transform of Eq. ([2|) which is done using the general rule 
for Fourier transforms 

F(x, t)^j^jd'kj doj e''^— *^^(k, uj). (4) 
In this case one finds 1191 



gL(k,cj) 
gT(k,w) 
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(5) 
(6) 
(7) 



where the source and perturbed momentum density vec- 
tors are divided into transverse and longitudinal parts: 
g — kg^ + gT and J = kJ^ + Jt, with k being the unit 
vector in the direction of k. Eq. (I7|) is a diffusion equa- 
tion and does not describe sound propagation. The quan- 
tity gT gives the diffusive momentum density generated 
by the fast parton. Eqs. ([5]) and ([6|), on the other hand, 
are damped wave equations. They describe damped sound 
waves propagating at speed Cs- Conical Mach cone struc- 
ture in the medium dynamics are generated by de and gi . 

The medium's dynamics can now be solved for by in- 
serting the explicit form of J'^ in Eqs. ©-([Tl) and Fourier 
transforming back to position space. As mentioned above, 
I will work with the source term derived in [TB] , where the 
fast parton was treated as the source of an external color 
field interacting with a perturbative QGP in the context of 
kinetic theory. In the relativistic limit (7 = (1 — u^)^^/^ ^ 
1), the position space representation of the source term in 
a gluonic medium at temperature T in the presence of a 
parton moving with velocity u = ux at position r — utx 
is given by 



rix) = {J"{x),uj"{x) 



where 



J°(p,z,t) = d(p,a:;,i)7u2 1^1- 
Jv(p, X, i) = (x - ut) d{p, x, t) 



and 



d{p,x,t) 



8^(p2+ 723.^)3/2- 



(8) 

(9) 
(10) 

(11) 



In the above expressions, p ~ [z"^ + V^Y^"^ is the radius 
transverse to the x axis, (Qp)^ = 3 if the source parton is 
a gluon and 4/3 for a quark, as = 5^/47r is the strong cou- 
pling, = (x — u£), and ttid = gT is equal to the Debye 
mass for gluons in the hard thermal loop approximation. 

The solutions for (5e(x, t) and g(x, €) in the presence of 
the source term, ([8]), for a wide range of parameters have 
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been presented in [17], and in greater detail in [19]. The 
resulting dynamics have been shown to contain a prop- 
agating Mach cone and diffusive wake. My purpose here 
is not to revisit those results, but rather to use them to 
construct a distribution function, which can then be used 
in the CF freeze-out prescription to obtain an azimuthal 
final particle spectrum. 

Given solutions for 6e(x,t) and g(x, t), one can de- 
termine the induced flow velocity and temperature mod- 
ification in the medium. The ideal hydrodynamic energy 
momentum tensor is given by 



(12) 



where u'^ is the flow velocity of the medium, which for 
the unperturbed static medium is just given by u'* = 
(1, 0, 0, 0). Thus, working to first order in perturbed quan- 
tities, one has for T^' 



g 



(eo + Se + po + Sp){ua + 6u)°{ua + SuY (13) 



Cooper Frye Brick Freezeout Scenario 




I conclude that 



Fig. 1. (Color online) The azimuthal particle spectrum gen- 
erated by a fast parton in an isochronous freeze-out scenario 
(as discussed in the text). The peak centered at (j> ~ has 
been added to simulate the near side distribution observed at 
RHIC, while the away side peak centered ai (j) — n results from 
the medium dynamics. Although the medium response consid- 
ered here is purely conical in structure, the away side peak is 
a broad Gaussian rather than double peaked. The isochronous 
CF freeze-out scenario has washed out any evidence of the con- 
ical structure present in the medium's dynamics. 



6u 



eo +Pq 



3g 

4eo 



(14) 



where I have used eo = 3po- One can express Su'^ in terms 
of 5u by requiring the medium's flow velocity remains 
physical, or u'^u^ = 1, which eventually yields 



Suq 



u • Su 

"00 



(15) 



An expression for ST can similarly be found by using the 
equation of state for an ideal gas of gluons 



15 



in which case one obtains 



4eo 



(16) 



(17) 



3 Freeze-out 

Having expressions for the flow velocity and temperature, 
it is now possible to construct the medium's distribution 
function, which in the Boltzmann limit, and up to nor- 
malization, is given by 



/(r,p,t) = e 



exp 



{up + Su^Pfj, 
To + ST 



(18) 



As discussed above, the distribution is converted into an 
azimuthal particle distribution by using a Cooper-Frye 



freeze-out scenario |23j . Consistent with the approach dis- 
cussed in P^ll^ , the final azimuthal particle spectrum for 
massless particles at mid-rapidity (y = 0) is given by 



dN , 



r dpTPT f dU^P^{f{r,p,t) - fo) 

(19) 



where Z"^ is the freeze-out hypersurface and 

= {pt,Pt cos(7r — (j)),pT sin(7r — <p),0). 



(20) 



The isotropic background contribution, /g = e~*"o , is 
subtracted in In what follows, it is understood that 
the source parton propagates along the —x axis. 

I now consider the spectrum given by (|19p . using the 
results for 5e{x,t) and g(x, t) which are discussed in |17j . 
The following parameters are used: = I/tt, Cs = 0.57, 
r]/s = l/47r, 7 = 33, (Q^)^ = 3 and T = 350 MeV. For 
the sake of clarity, in what follows I will not include the 
diffusive contribution to g from gx- This will ensure that 
all contributions to the final particle spectrum come from 
the conical Mach cone structure generated by Se and g^ . I 
have also added a Gaussian ansatz to the distribution cen- 
tered around (/> = to mimic the near side, or trigger, con- 
tribution to the hadron correlation functions observed at 
RHIC. The physical contribution, that obtained from the 
solutions for (5e(x, t) and gL(x, t), will be centered around 

(f) = TT. 

I first consider an isochronous freeze-out scenario, sim- 
ilar to that done in |22| . In this scenario, the hypersurface 
is given by dSf^ = (dF, 0,0, 0). The matter is assumed 
to instantly convert into free streaming particles, and the 
hypersurface integration in (|19p becomes an integration 
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Cooper Frye Radial Freezeout Scenario 
Radial Flow = 0.6 c 
30| ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 ■ r 




Fig. 2. (Color online) The azimuthal particle spectrum gen- 
erated by a fast parton in a radial shell freeze-out scenario 
with artificially added background velocity (as discussed in the 
text). The away side spectrum contains a markedly different 
structure than that seen in Fig. [T] In particular, the away side 
distribution now contains a double peak structure, reflecting 
the conical flow present in the medium. 

over all space. For the momentum integration in I 
have chosen = 4 GeV and ptp — GeV. The result is 
shown in Fig. [T] for an arbitrary normaUzation. As men- 
tioned above, the contribution around = simulates the 
near side distribution and has been added to provide con- 
text. The away side distribution is centered around (f) ~ n 
and is generated by the solution for (5e(x, t) and gL(x, t). 
Even though Se{x,t) and gL(x, i) only contribute a coni- 
cal Mach cone structure to the medium's dynamics, there 
is no sign of this in the away side distribution, which has 
a broad Gaussian shape instead of a double peak struc- 
ture. The isochronous CF freeze-out scenario has washed 
out any evidence of the conical structure present in the 
medium's dynamics. 

I next consider a scenario in which the freeze-out sur- 
face is a radial shell extending from /? = 4 — 5 fm. The 
source parton is assumed to start just inside the shell and 
then propagate through it until any contribution from the 
medium's hydrodynamic response is negligible. Recall that 
p = {z^+y^y^"^ and the source parton is propagating along 
the —X axis. The hypersurface is now oriented in the pos- 
itive p direction and the integration includes time. I limit 
the integration to a thin slice in the 2 direction, extend- 
ing from ±0.1 fm. Finally, an artificial radial background 
flow velocity of Vp — 0.6 is introduced in the medium ve- 
locity appearing in p9p . The result is presented in Fig. 
[21 where again, the normalization is arbitrary. As one can 
see, under this freeze-out scenario the away side distribu- 
tion now contains a double peak structure, reflecting the 
conical flow present in the medium. 



4 Discussion and Conclusions 

The results presented above show that theoretical cal- 
culations of the azimuthal away side hadron correlation 



spectrum generated by fast partons are sensitive to the 
specific freeze-out scenario employed. Specifically, it was 
shown that an isochronous Cooper-Frye freeze-out sce- 
nario yielded significantly different results than a radial 
shell scenario with an artificially added background veloc- 
ity (as discussed in Sec. [3]). The results, as plotted in Figs. 
[1] and [21 are not meant to provide a meaningful compari- 
son with the data observed at RHIC, but are rather shown 
to highlight the sensitivity of results to the freeze-out sce- 
nario chosen. 

Indeed, even the result presented in Fig. [2] is sensi- 
tive to the choice of parameters. A change in the speed 
of sound in the medium from Cg = 0.57 to Cg = 0.3 has 
the effect of wiping out the double peak structure, leaving 
a single Gaussian peak on the away side. A similar effect 
is observed if one changes rj/s from l/47r to 6/47r while 
leaving the speed of sound unchanged. 

The potential for ambiguity in various freeze-out sce- 
narios is not surprising. Freeze-out occurs in physical sys- 
tems when the mean free path goes to infinity and particles 
become free-streaming. This occurs in heavy-ion collisions 
because the medium is expanding and cooling. A conven- 
tional choice for the freeze-out temperature is 110 MeV 
|25j . well below the initial temperature of the medium 
(around 400 MeV) and the critical temperature for hadroniza- 
tion (around 160 MeV). An infinite static background 
medium, such as that considered above and in the compar- 
ison done in [22], does not expand or cool, hence there is no 
reason for freeze-out to occur. Freeze-out must be imposed 
on the system independent of physical processes, creat- 
ing the potential for widely varying and inaccurate re- 
sults. Any meaningful theoretical comparison of azimuthal 
hadron correlation functions measured at RHIC requires 
implementing a realistic freeze-out scenario in an expand- 
ing medium. 
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